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The  analysis  of  wave  propagation  in  periodic  elastic  structures, 
i.e.,  structures  in  which  the  elastic  constants  of  the  material  vary 
• in  a periodic  manner,  is  of  great  interest  due  to  the  presence  of  passing 

and  stopping  frequency  bands  in  such  structures.  While  such  an  analysis 

if  * 

is  in  general  quite  complicated,  closed  form  solutions  have  recently 

■ 

been  obtained  for  a number  of  simpler  problems.  It  is  the  purpose  of  this 
note  to  discuss  an  alternative  approach  to  such  problems  which  may  be  of 
some  value  in  gaining  physical  insight  into  the  properties  of  the  solu- 
tion and  in  qualitatively  "sketching"  the  frequency  spectrum  with  reduced 
effort. 

One  common  physical  case  studied  so  far  is  that  of  torsional 
elastic  wave  propagation  in  a bi-element  composite  circular  cylinder, 
either  solid,  reference  [1],  or  hollow,  reference  [2].  The  outer  (and 
inner,  if  present)  boundaries  are  taken  as  stress  free.  The  prismatic 
cylinder  is  assumed  to  be  infinite  in  length,  |z|  < <*>,  and  consists  of 
two  homogeneous  elements,  comprising  region  I of  length  £,  density  p 
and  shear  modulus  p and  region  II  of  length  density  p'  and  shear 
modulus  p',  repeated  periodically  in  the  axial  (z)  direction  with  a unit 
cell  size  d = (Jt  + 4').  Such  a problem  has  only  a single  displacement  com- 
ponent, U , independent  of  the  0 coordinate  in  a cylindrical  coordinate 

system  (r,0,z).  This  displacement  component  satisfies  a wave  equation 

1/2  1/2 

with  phase  speed  c = (p/p)  and  c'  = (p'/p')  in  regions  I and  II 
respectively.  We  choose  the  z coordinate  positive  to  the  right  and  assume 
a time  harmonic  dependence  proportional  to  exp(-iwt),  where  u is  the  ang- 
ular frequency  in  radians  per  unit  of  time. 

In  this  simple  example,  the  governing  field  equation  readily  separates 
in  (r,z)  to  give  tangential  component  of  displacement  in  the  form: 


(1)  (a)  Region  I,  0 < z < 4: 

■ [A  sin  C Az/4)  + B cos  (Az/4)]  Z^Cxr/a)  exp(-iwt), 

0 

Cb)  Region  II,  - 4*  < z < 0: 

uP  = [A'  sin  (A'z/4')  + B cos  (A’z/4')]  Z^Pic'r/a)  exp(-iwt). 

We  take  the  outer  radius  of  the  cylinder  to  be  a,  the  inner  radius 

(if  one  exists)  b,  the  half  thickness  h = (a-b)/2  and  ft  = w/ojs  where 

<*>  = 7rc/2h  is  the  lowest  thickness  shear  frequency  of  an  infinite,  iso- 

s 

tropic,  homogeneous  plate  of  half  thickness,  h,  mass  density,  p,  and  shear 
modulus,  y.  We  have  condensed  the  notation  (which  corresponds  to  [1] 
and  [2])  by  introducing  k e irka/2h,  k'  = irk*a/2h,  X = ir£i/2h,  X'  = ir5'4'/2h. 
The  radial  solutions,  Z,  are  given  by 

(2)  Z1  « (2/k)  J1(lcr/a)  <2  >_  0 

■ (2/k)  I^(tcr/a)  <2  <_  0 

- r/a  k2  * 0 

for  the  solid  cylinder  and 

(3)  Z1  - (2/k)  J^Kr/a)  - (wk/2)  DY^Kr/a)  <2  > 0 

■ (2/k)  I^(xr/a)  + k D K^(Kr/a)  k2  <_  0 

2 

■ r/a  + Da/r  < “0 

for  the  hollow  cylinder  with  similar  equations  for  in  terms  of  <' 


and  D'  and  where 
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(4)  Uic/a)2  + X2  - Uco/c)2  or  k2  + ^ - O2, 

U'K’/a)2  + X'2  - UVc')2  or  k'2  + £'2  - (c/c')2fi2. 

The  radial  boundary  conditions  require  zero  shear  stress,  x =» 

i 

Vir3[U0/r]/3r,  on  the  outer  surface  r=a  and  either  zero  shear  stress  on 
the  inner  surface  r=b  for  the  hollow  cylinder  or  a bounded  displacement 
at  the  origin  0 for  the  solid  cylinder  (this  was  already  implied  in 
the  form  of  Z used  above  for  the  solid  cylinder) . 

For  the  solid  cylinder,  the  radial  boundary  condition  requires 

(5)  J2(k)  = 0, 

2 2 

for  ic  >_  0 and  has  no  roots  for  k < 0.  Roots  of  this  transcendental 
equation  as  given  in  [3]  are  Kq  * 0,  vc^  - 5.1356,  <2  * 8.4172, 
etc.  »- 

For  the  hollow  cylinder,  the  boundary  conditions  require 

(6)  J2(k)  Y20ct)  - J2(xt)  Y2(k)  =0,  t = b/a 

D = (ttk2/4)  Y2(k)/J2(k), 

2 2 
for  ic  >0  and  no  solution  for  k <0.  Roots  of  this  transcendental 

equation  must  be  calculated  as  a function  of  t.  For  t ■»  1/3,  - 

1.1892,  «2  = 2.1160,  = 3.0811,  etc.,  and  for  t = 1/2,  *=  1.08446, 

<2  ■ 2.04602,  ■ 3.03122,  etc.  Similar  results  hold  for  k*  . 

The  remaining  boundary  conditions  apply  to  surfaces  of  constant  z. 

The  origin  and  direction  of  the  z axis  are  arbitrary;  we  choose  z-0 

at  the  interface  between  regions  II  and  I and  first  require  continuity  of 

Ug  and  t ■ p3Ug/3z.  This  leads  to  the  two  conditions 


i|{> 
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(7)  vA  - v'A',  B - B', 

where 

v = yA/A.  and  v'  = y'A'/A  . 

The  remaining  interfacial  condition  is  a quasi-periodicity  condition 
which,  by  Floquet's  theorem,  requires  the  wave  amplitude  to  have  the  same 
periodic  structure  as  that  of  the  elastic  medium,  in  order  to  have 
a solution 

U0(r,z,t)  = w(r ,z)  exp  1 (yz  - u>t), 

where  w is  periodic  in  z with  periodicity  d,  i.e.,  w(r,z+d)  » w(r,z) 
and  y is  the  Floquet's  wave  number  corresponding  to  the  phase  shift, 
which  has  to  be  determined  from  the  solution  of  the  problem. 

This  leads  to  the  quasi-periodic  conditions 

(a)  U0  (r.A.t)  = Ug1  (r.-A'.t)  exp  iyd, 

\8) 

(b)  t^0  (r,A,t)  = (r,-A',t)  exp  iyd, 

where  for  convenience  we  have  dropped  explicit  time  harmonic  dependence 
exp(-iut).  Equation  (8)  implies  Z^icr/a)  = ZII(ic'r/a)>  i.e.,  k = k ' ; 
and 

(a)  A sin  \ + B cos  A ■ [-A'  sin  A'  + B'  cos  A']  exp  iyd, 

(9) 


(b)  vlA  cos  A - B sin  AJ  * v'[A'  cos  A'  + B'  sin  A']  exp  iyd. 
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While  the  problem  is  completely  formulated  at  this  point  and  has 

indeed  been  solved  in  this  form  in  references  [1]  and  12],  we,  however, 

introduce  slight  modification  in  the  variables  of  the  problem. 

We  introduce  the  "interfacial  elastic  parameter",  Q,  as  the  ratio 

of  stress  to  displacement  at  an  interface.  Clearly  Q is  piecewise 

continuous  throughout  the  rod.  Consider  first  q|  = Q defined  as 

ZBU  o 


(10)  Qo  - Tz0  (O)/U0(O)  - VA/B  - v'A'/B'. 


Equations  (9a)  and  (9b)  may  then  be  rewritten  as 


(a)  B(Qq  sin  X + v cos  X)  = B'(v/v')  [-Q  sin  X'  + v'  cos  X']  exp(iyd). 


(b)  B(Qo  cos  X - v sin  y ) = B ' [ Qq  cos  X'  + v'  sin  X']  exp(iyd). 


Consistency  of  these  two  homogeneous  equations  requires  that  the  deter- 
minant of  the  coefficients  be  zero.  This  leads  us  to  the  frequency 
equation 


2 2 2 
Q [v'  sin  X cos  X'  + v sin  X'  cos  X]  + Q (v'  - v ) sin  X sin  X' 

o o 


(12)  + vv'  (v'  sin  X'  cos  X + v sin  X cos  X')  - 0,  0<Q  <°°. 

— o 


On  the  other  hand,  we  may  eliminate  Qq  from  Eq.  (11)  by  using  Eq.  (10) 
and  then  obtain  the  usual  dispersion  equation. 


[exp(iyd)  ] + [(v/v'  v v'/v)  sin  X sin  X'  - 2 cos  X cos  X'J  exp(iyd)  +1-0, 
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which  relates  Floquet  wave  number  y to  the  axial  wave  numbers  X and  X'. 
Clearly,  the  second  order  equation  has  two  solutions  whose  product  Is 
unity  and  whose  sum  is 


(14)  cos  yd  **  cos  * cos  X'  - — (~t  + ) sin  X sin  X', 


If  we  examine  the  completely  homogeneous  case,  H = p = p', 
p ■ p',  £ = £'  and  v = v',  we  get  the  solutions  y *=  ± (ir/d)  (££/h  ± 2n), 
n = 0,1,2,...  . On  the  extended  zone  scheme  we  select  n = 0,  and  therefore 
the  two  values  for  y = ± 7r£/2h  = ± X/S.  are  real  and  simply  represent 
the  two  directions  in  which  torsional  waves  can  propagate  with  no  change 
in  amplitude.  This  is  best  seen  by  examining  the  form  of  the  z dependence 
in  the  tangential  component  of  displacement 


U.  * B[ (A/B)  sin(Xz/£)  + cos(Az/£)]  Z(<r/a), 
0 


where,  from  equation  (10)  it  can  easily  be  shown  that 


A/B  = Qq/v  = ±i. 


Hence,  as  expected 

U0  - B[cos (Xz/i)  ± i sin(XzM)]  Z(tcr/a), 


■ B exp(±iXz/i)  Z(xr/a), 
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Whtle  there  may  be  no  general  advantage  In  calculating  Qq  for  a given 
fJ  and  then  calculating  y etc.,  there  are  some  special  circumstances  where 
the  Introduction  of  Qq  is  useful.  Consider  first  the  limiting  cases  when 
Qq  tends  to  0 or  <*>,  corresponding  to  zero  stress  or  zero  displacement, 
respectively,  at  z=0.  The  former,  Qq=0  leads  to  (if  cos  X and  cos  X' 


are  non- 

zero) 

Ca) 

v tan  X 

+ V 

' tan  X'  = 0,  v,v' 

' ^ 0 

(b) 

cos  yd  1 

1 

B — 
2 

(cos  X'/cos  X + cos 

X/cos  X')» 

(15) 

(c) 

uo  ■ B ' 

cos 

(Xz/i.)  Z(icr/a) , 

(d) 

““  - B 

cos 

(X'z /l*)  Z(Kr/a). 

Similarly,  for  Qq  = 00 , we  have  (if  sin  X and  sin  X'  are  non-zero) 

(a)  v cot  X + v'  cot  X'  = 0,  v,v*  / 0 

(b)  cos  yd  = y (cos  I'/cos  X + cos  X/cos  X'), 

(16) 

(c)  Ug  - A sin  (Xz Jl)  Z(icr/a), 

(d)  U*1  = A(v/v * ) sinCX'z/i')  Z(icr/a). 

0 

These  frequency  equations  for  Qq  ■ 0 and  ®,  are  not  only  simpler  to  solve 
than  the  general  set,  they  are  formally  of  the  same  structure  as  the 
frequency  equations  for  the  endpoints  of  the  Brillouin  zones  which  define 
the  passing  and  stopping  bands.  Equations  at  the  zone  ends  are  obtained 
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by  setting  exp(iyd)  - ± 1 in  equation  (13);  these  are  the  cut-off  points 
where  y changes  from  real  to  Imaginary  values.  For  expClyd)  ■ + 1» 
this  has  solutions  with  a period  d at  frequencies  defined  by 

(a)  v tan(X/2)  + v'  tan(X'/2)  - 0, 

(17) 

(b)  v tan(X'/2)  + v'  tan(X/2)  - 0, 

and  for  exp(iyd)  = -1,  it  has  solutions  of  period  2d  at  frequencies 
defined  by 

(a)  v tan(X/2)  - v'  cot(X'/2)  = 0, 

'b)  v cot(X/2)  - v'  tan(X'/2)  = 0. 

We  thus  see  that  the  wave  numbers  satisfying  equation  (15(a))  are 
one-half  that  satisfying  equation  (17(a));  that  for  equation  (16(a))  are 
one-half  the  wave  numbers  satisfying  equation  (17(b)). 

It  may  be  pointed  out  that  satisfies  a quadratic  equation,  and 
for  a fixed  ic  and  given  values  of  (S},y)>  there  are  always  two  values  of  Qq. 
Thus 

(18a)  Q^1}  = 0;  = v tan  X = -v'  tan  X', 

are  the  two  roots  of  Eq.  (12)  in  the  presence  of  the  constraining  equa- 
tion (15a) . The  displacement  field  in  the  second  case  is 

Ug  » A c*c  X cos  X(z/Z  - 1)  Z(icr/a), 

(19a) 

Ug1  ■ -A'  esc  X’  cos  X * (z/ S. ’ + 1)  Z(icr/a). 
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Stmilarly 


(18b) 


.CD 


» ^o^  - -v  cot  A » v'  cot  A', 


are  the  two  roots  of  the  same  equation  in  the  presence  of  the  constraining 
equation  (16a).  The  displacement  field  in  the  second  case  is 


(19b) 


Uq  ■ A sec  A sin  A(z/£  - 1)  Z(icr/a), 


Ug1  - A'  sec  A'  sin  \'(z/V  + 1)  Z(»cr/a). 


In  each  of  these  four  cases  the  Floquet  number  y is  governed  by  the  same 
equation  (15a). 

These  "second"  solutions  correspond  to  the  same  problem  with  the  z 
origin  shifted  by  one  segment,  e.g.  by  £,.  For  example,  the  solution  for 
“ 0 corresponds  to  a zero  shear  stress  at  z equal  to  zero  while  the 


(2) 


"second"  solution,  Qq  = v tan  A,  corresponds  to  a zero  shear  stress 
at  z = £,,  etc. 
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The  equation  (15b  - 16b)  can  be  rewritten  as 


(20) 


sin(yd/2)  = ± ((cos  A/cos  A')^^  - (cos  A'/cos  A)^^), 


cos(yd/2)  = ± y ((cos  A/cos  A')*^  + (cos  A'/cos  A)^^), 


and  therefore  we  see  that  at  the  end  points  of  the  Brillouin  zones 


(21) 


cos  A'  - cos  A = 0 when  yd  = 0,  2ir,  4tt, 


cos  A'  + cos  A = 0 when  yd  = ir,  3ir,  5tt, 


Thus  on  the  left  end  of  the  zone 


(22a)  A'  = A ± ne  ° 0,  2,  4,  ... 


and  on  the  right  end  of  the  zone 


(22b)  A'  =>  A ± noff,  nQ  = 1,  3,  5,  ... 


Consider  now  the  case  when  Q =0.  On  the  left  end  of  the  zone 

o 

cos  A'  = cos  A imply  sin  A'  = sin  A and  therefore  from  Eq.  (15a) 


(23)  (v  + v')  tan  A = 0. 


Consequently,  either  (v  + v')  = 0,  or  sin  A = 0.  Hence,  either 


(24)  A - inn^/Kp/w)  (I/I)  + 1],  or  A = mir,  m = 0,  1,  2, 


Knowing  A from  this  equation,  we  can  now  find  A'  from  Eq.  (22a).  However, 
in  the  second  case,  n in  Eq.  (22a)  and  m.  in  Eq.  (24).,,  cannot  both  be 
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arbltrary  since  from  Eq.  (4) 

(25)  X'2  » U'/i)2 (c/c)2  X2  + U’/a)2  2 CCc/c')2  - 1). 

When  = °°,  we  find  from  Eqs.  (21 ) ^ and  (16a)  that 

(v  + v')  cot  X = 0, 

and  therefore,  either  (v  + v')  = 0,  or  cos  X = 0.  Hence,  in  this  case 
either 

(26)  X = ±Trne/[(p/p)(i’/i)  + 1],  or  X = mir/2,  m = 1,  3,  5,  ... 

On  the  right  end  of  the  zone,  we  have  simply  to  replace  the  integers 

by  rio  to  obtain  corresponding  formulas  for  X and  X'  when  Qq  = 0 or 

Q *=  00 . 

o 

It  is  easy  to  see  that  in  the  interior  of  the  Brillouin  zone  0 < yd  < it, 
yd  is  always  complex  for  real  values  of  X and  X',  X ^ X'.  This  follows 
immediately  from  the  fact  that  Eq.  (15b)  can  be  rewritten  as 

(27)  exp  iyd  = cos  X/cos  X'. 

We  therefore  conclude  that  when  Q = 0 (stress-free  interface) , or  when 

c 

Q ■=  00  (displacement-free  interface)  torsional  waves  cannot  be  propagated 
o 

in  a prismatic  cylinder  with  periodic  structure  and  are  therefore  damped 
out. 

We  now  consider  those  values  of  Q whifh  lead  to  information  concern- 

o 

ing  frequencies  in  the  passing  band.  From  Eq.  (12)  we  find  that  when 

Q - ±iv' 
o 


For  every  value  of  radial  wave  number  k,  these  two  equations  determine  a 

unique  value  of  ft  and  y.  Plotting  ft  versus  yd,  the  first  equation  gives  us 

lines  of  constant  ft  for  every  given  k and  for  different  values  of  n- 

The  second  equation  gives  us  a series  of  hyperbolas  (or  straight  lines  when 

k^  = 0).  Starting  with  the  cut-off  frequencies  when  yd  = 0,  the  real 

branches  of  the  dispersion  curves  lie  between  the  bounds  defined  by 

Eq.  (31).  The  dispersion  curve  will  cross  the  bounds  only  at  their  points 

of  intersection,  where  Q = ± iv'. 

o 

Another  independent  set  of  bounds  can  similarly  be  found  when 
Qq  = ± iv,  leading  to  sin  X’  =0.  In  this  case  the  two  equations  of  the 


bounds  are 
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(c/c')V  - k2  + C2nh/Z')2, 


fi2  - k2  + C2h/Trj02(Yd  ± mrr)2. 

Since  the  cut-off  frequencies  at  y = 0 and  yd  = it  are  available  from 
the  roots  of  Eqs.  (17,  18),  the  knowledge  of  additional  intersection  points 
inside  the  Brillouin  zone,  provides  us  with  sufficient  information  to 
sketch  the  dispersion  curve  qualitatively.  The  use  of  interfacia!  parameter 
Qq  therefore  provides  us  with  a scheme  for  a qualitative  solution,  which 
may  be  valuable  in  more  complex  problems 

Physically,  these  two  cases  correspond  to  A’  = iB'  and  A = iB 
respectively.  These  conditions  in  turn  lead  to  the  requirement  that  the 
ratio  of  shear  stress,  tz0(z),  to  displacement,  UQ(z),  be  uniform 
throughout  region  II  for  Qq  = iv  and  throughout  region  I for  Qq  = iv', 
i-e*»  t20II(z^/U0II(z^  = ^g^COJ/U^CO)  respectively . 

The  bounding  curves  for  Qq  = iv,  and  iv',  are  shown  in  Fig.  1 along 

with  the  actual  dispersion  curve  as  calculated  in  ref.  [1]  for  the  lowest 

branch,  = O.for  a solid  cylinder.  The  parameters  used  are  p/p’  = 1/40, 

l/a  = 3/5,  Jl'/a  = 3,  c/c'  = 1/4.  For  these  parameters,  the  bounding 

vd  3 

curves  for  Q = iv'  are  0 = 5n/3:  -1 — =7-0  + 11  and  for  Q = iv  are 
o rr  4 o 

0 = 4n/3;  a + n. 

TT  J 

The  scale  for  0 in  reference  [1]  is  different  from  that  used  in 


this  paper  by  a factor  of  2 since  u>  there  was  nondimensionalized  with 
respect  to  the  lowest  thickness  shear  mode  frequency  of  an  infinite  plate 
of  half  width  a while  here  a half  width  h,  equal  to  a/2  for  a solid 
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El 


cylinder,  is  used.  Thus  values  from  reference  [1]  used  for  comparison 

here  must  be  divided  by  2.  Thus,  the  cut  off  frequencies 

are  at  !5  ■=  0,  0.290,  1.168,  1.359,  1.628,  1.872,  2.620,  2.772,  etc. 

The  intersection  points  of  the  bounding  curves  lie  for  Qq  = iv  at  (1.8, 

1.33)  for  n = 1,  m = 1 and  at  (3.6,  2.67)  for  n = 2,  m = 2 while  for 

Qq  = iv',  the  first  point  lies  at  (2.25,  1.67)  for  n = 1,  m = 1 and  the 

others  lie  beyond  the  range  of  calculated  values.  The  complex  roots 

lie  at  (1  + 0.177i,  0.764)  and  (1  + 0.729i,  0.679)  for  Q = 0,  °° 

o 

respectively  on  the  first  complex  branch,  at  (2  + 0.055i,  1.361)  and 

(2  + 0.080i,  1.1624)  for  Qq  = 0,  " respectively  on  the  second  branch,  etc. 

The  effort  in  calculating  all  of  these  bounds  and  solution  points 
and  the  end  points  of  the  Brillouin  zones  is  considerably  less  than  that 
involved  in  determining  the  complete  dispersion  spectrum,  even  for  this 
simple  problem,  yet  these  few  calculations  are  sufficient  for  a reasonably 
accurate  sketch.  It  is  anticipated  that  this  reduction  in  effort  will 
apply  to  more  general  problems  as  well. 
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